Clones and other interference effects in the evolution of angular momentum coherent 

states 



The aim of this article is to present the interference effects which occur during the time evolution 
of simple angular wave packets (WP) which can be associated to a diatomic rigid molecule (het- 
eronuclear) or to a quantum rigid body with axial symmetry like a molecule or a nucleus. The time 
evolution is understood entirely within the frame of fractional revivals discovered by Averbukh and 
Perelman since the energy spectrum is exactly quadratic. Our objectives are to study how these 
interference effects differ when there is a change of the initial WP. For this purpose we introduce 
a two parameter set of angular momentum coherent states. From one hand this set emerge quite 
naturally from the three dimensional coherent states of the harmonic oscillator, from another hand 
this set is shown to be buit from intelligent spin states. By varying one parameter (77) a scenario 
of interferences occur on the sphere at fractional times of the revival time that strongly depends 
on 77. For 77 = ±1 the WP, which coincides with a WP found by Mostowski, is a superposition of 
Bloch or Radcliffe's states and clone exactly in time according to a scenario found for the infinite 
square well in one dimension and also for the two dimensional rotor. In the context of intelligent 
spin states it is natural to study also the evolution by changing 77. For 77 = the WP is called linear 
and produces in time a set of rings with axial symmetry over the sphere. The WP for other values 
of 77 are called elliptic and sets of fractional waves are generated which make a transition between 
two symm etries. We call 'mutants' these fractional waves. For specific times a clone is produced 
that stands among the mutants. Therefore the change in 77 produces novel change in the quantum 
spread on the sphere. We have also constructed simple coherent states for a symmetric rotor which 
are applicable to molecules and nuclei. Their time evolution also shows a cloning mechanism for 
the rational ratio of moments of inertia. For irrational values of this ratio, the scenario of partial 
revivals completed by Bluhm, Kostelecky and Tudose is valid. 



In the recent years one has found the existence of a generic behaviour for the time evolution of simple quantum 
systems. In ref. (I]] Averbukh and and Perelman have indeed discovered a universal scenario of fractional revivals in 
the long term evolution of quantum wave packets of bounded system which goes beyond the correspondence principle. 
They established this scenario by expanding the bound states energies relevant to the wave packet up to the second 
order with respect to the mean energy, thus producing a local spectrum linear plus quadratic in one quantum number. 
For a well concentrated wave packet at the initial time they defined two time constants T c i and T rev > T c i such that 
the foregoing evolution of the wave packet is predicted as follows: for < t < T c i, the wave packet spreads around 
the mean trajectory that can be associated with the underlying classical evolution while for T c ; < t < T rev the wave 
packet interferes with itself in such a way that for fractional times t — (m/ 71) T rev the wave packet is divided into 
q fractional wave packets. If n is even, q = n/2, if not, then q — n. The specificity of the system and that of the 
wave packet determines the shape of the fractional wave packets which are supposed to be spread regularly around 
the mean trajectory. At time t — T rev , the wave packet is rebuilt either identical or only resemblant to the original 
one, depending on the importance of the neglected terms higher than quadratic. In some very specific cases, the 
fractional wave packets are clones of the initial one as studied in the most recent paper Q and the revival is exact. 
This scenario has been validated in the most spectacular manner by several authors for a wave packet in a circular 
orbit of the hydrogen atom ||, see also [Q and by Boris et al. || for an elliptic orbit. It has been extended by 
Bluhm and Kostelecky Q to a very long evolution with the demonstration of superrevivals due to the cubic terms. 
The analytical explanation of the effects caused by cubic terms has been found by Leichtle et.al. j7j. Examples of 
vibrational wave packets in anharmonic potential of simple molecule were also found since [|J,[|. The scenario has 
been thereafter extended to cases where the energy depends on two quantum numbers p0| . A recent synthesis |Tl] ] 
contains most of the references on this topic while ref. jFq] has been devoted to the definition of the experimental of 




Theoretical Physics Department, Institute of Physics, 
University Maria Curie-Sklodowska, PL 20-031 Lublin, Poland 

and 

Institut des Sciences Nucleaires, F 38026 Grenoble- Cedex, France 



R. ArvietE 



Institut des Sciences Nucleaires, F 38026 Grenoble-Cedex, France 
(Received 13 January 1998, in revised form 10 July 1998) 



I. INTRODUCTION 
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wave packets in atomic physics and in molecular physics. The recurrences studied in ref. |1|-|12|] are all devoted to wave 
packets and can be mathematically explained in terms of special Gauss'sums. It is necessary to point out that such 
sums were earlier extensively used in the physical litterature by Berry and Golberg |p|, and Berry Jl3] who studied 
the full time evolution of the propagator of a nuclear spin with an hamiltonian of the form l 2 /2J. These authors 
developed a renormalization theory and discussed the semiclassical limit. More recently M| various Talbot effects 
(integer, fractional, fractal) were also discovered with the help of these techniques. Finally Berry has shown the 
occurrence of fractal dimensions both in space and time during the evolution of a uniform WP in boxes of arbitrary 
dimensions. 

Our aim is to consider only angular wave packets for diatomic molecules and also for symmetric rotors and our 
purpose is to discuss the time evolution of a large enough set of coherent states of the angular momentum. There 
is rather extensive literature on the coherent states. A full bibliography on this subject is found in r cf. Q. But 
we will focus on the intelligent spin states described in ref. |18|-|23| Large efforts have been made to build up such 
wave packets mainly in the context of group theory. Little effort, on the contrary, has been made to understand in 
detail their time evolution. Until the work of jy it was not realized that they could evolve according to an universal 
scenario. For a diatomic molecule or for a symmetric rotor, the spectrum of which are quadratic in quantum number 
like for the systems discussed in ^ , the universal scenario is exact and repeats with period T rev . The main 
question is whether the fractional wave packets are clones or just resemblant to the initial ones. Despite all efforts 
to concentrate initially the wave packet in the best possible way, the quantum evolution destroys this concentration 
(spatial localization) according to the rule formulated in |l[ . We will show in this paper under which conditions the 
wave packets separate into clones and in which conditions there is a more restricted scenario of partial revivals. We 
will also show the existence of new fractional waves with different shapes that we call mutants. It is crucial to present 
a simple and physically meaningful picture of a coherent state and to decrease the number of its parameters to the 
minimum. One should keep in mind that a coherent wave packet has a classical content larger than the eigenstates 
of the angular momentum. Due to this property, it is possible to choose in the simplest manner the coherent wave 
packet as we will show below. 

In section || we will study a large set of coherent angular WP which depend only on the angles 8 and <f>. If we 
impose one condition of minimum uncertainty these WP are composed of eigenstates of L 2 and contain, after a proper 
choice of the axis of coordinates, only spherical harmonics with magnetic quantum number m of the same parity as I, 
These states belong to the family of coherent states called the intelligent spin states p^-|2^]. The main body of this 
article is organized around a particular subset: the exponential WP which are shown to be narrowly related to the 
coherent states of the harmonic oscillator (see the Appendix). Moreover the angular spread of the probability density 
depends on a single adjustable parameter. One of the limiting case is a WP derived originally by Mostowski |25| for 
a diatomic molecule and is called circular. By varying one parameter one obtains new WP with cylindrical symmetry 
that we call linear and a large set called elliptic. Each element of this set corresponds to a quantum system like a 
diatomic molecule or a nucleus in a pure state with a specific preparation giving to it an average angular momentum, 
an angular distribution and particular spread of the distribution of the angular momentum. 

In section [II it will be shown that a WP obtained by Atkins and Dobson |^(| using Schwinger j2?| boson repre- 
sentation of spin 1/2 is a particular circular state almost coincident with the exponential WP. More generally the 
boson representation for a spin s is also shown to lead to circular states. The exponential WP provides a closed and 
compact expression in the angular variables. 

In section |y| we will state precisely for our angular WP the scenario of fractional revivals derived in [[[] for the 
general case. It is indeed possible to specify the shape of the fractional revivals. The most spectacular event of 
cloning found essentially for the infinite square well in one dimension [^) and for a two dimensional rotor in |Tl| l 
is extended here to the most general case of circular WP. For the linear and elliptic WP the fractional waves are 
generally different from the initial WP due to the quantum spread. However for a particular set of times a single 
clone exists which coincides with the initial WP. In the case of the most general linear WP the fractional waves keep 
cylindrical symmetry showing isotropy in the spread over the sphere. For the elliptic WP the fractional waves should 
accomodate two limiting symmetries: the plane symmetry present for circular states and the cylindrical symmetry 
valid for the linear states. We call 'mutants' these intermediate fractional waves. 

In section we present a numerical calculation showing this evolution for exponential WP. The cases where cloning 
occurs is somewhat obvious, however it is interesting to define properly the time windows during which a given system 
of clones governs the time evolution. The 'carpet' representation used also elsewhere [f38| is an interesting tool in this 
respect. A second interesting result of this section lies in the shape of the mutants which can hardly be found from 
analytical considerations. It is found that these fractional waves preserve a good angular localization on the sphere. 
Their shape differs however from that of the initial WP since the generic structure is a well defined crescent-like shape. 

Finally in section VI , we will construct an angular momentum coherent state of a symmetric rigid rotor according 
to the rules defined by Janssen p8| . The time evolution of such a state is studied in section VII. Once the number 
of parameters is reduced, the time evolution of the coherent state, which is now a three dimensional system with two 
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quantum numbers as in ]iO| , presents clones if the ratio of the moments of inertia is rational. In the case when this 
ratio is not a rational number the fractional wave packets are not clones. 

II. DERIVATION OF COHERENT ANGULAR WAVE PACKETS 

Coherent angular WP can be defined as functions of 9 and <j> which fulfill two requirements: 

1. Their angular spread should be under control i.e. it should be possible to adjust their angular distribution in 
the easiest manner by changing a few parameters or a simple function. 

2. A criterion of minimum uncertainty should be obeyed in a manner similar to the conditions satisfied by the 
coherent states of the harmonic oscillator. 

AqiApi = ^, i = x,y,z . (1) 

In that direction many attempts have used uncertainty relations derived with the angular variables and the 
angular momentum operator ( see for example |fj"7f for a complete reference to these works). Despite the long 
list of works devoted to this field there is no detailed evolution of the time evolution of these WP in the literature 
as said our introduction. The first attempts in the line of the modern developments has been done in in 
which a WP of the rigid rotor in 2 dimensions was shown to clone exactly. Our work will extend these results 
to three dimensions and will discuss a rather large class of WP. 

In the following we will use only uncertainty relations based upon the components of the angular momentum and we 
will call coherent states all the states which satisfy 

ALlALl = ±{L z ) 2 . (2) 

We assume that this equation holds with the axis of coordinates such that 

{L x )=0, (L„)=Q. (3) 
It is a textbook result that a general WP satisfies the inequality (in units with % = 1) 

ALtALl>±\{[L x ,L v })\ 2 . (4) 

This result is derived by considering the norm of the state obtained by application of a special combination of L x 
with L y which involve a real parameter called rj 

(L x +ir,L y )\V) . (5) 
If the minimum uncertainty condition is realized there exists a value of rj for which 

(L x +i V L y )\^) = 0. (6) 
This value of r\ is related to the average values by two formulas 129] 



V 2ALI ± V AL y ■ (7) 

The second of these equations provides a meaningful interpretation of rj in terms of AL\ and ALy. Let us now 
construct states which satisfy eq. (||). 
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A. Eigenstates of L 2 . 



The simplest and most natural possibility is to construct the states \ty) as eigenstates of L 2 . This problem was 
solved long ago ]T^-|2^| and the solutions were called intelligent spin states. Let us briefly sketch a few of their 
properties and explain why we will consider only a subset of them. 

The intelligent spin states are the eigenstates \w) of L 2 and of the (non hermitian) operator L x + \r\L y with 
eigenvalue w such that 

(L x + vqLy) \w) = w \w) . (8) 



The (21 + 1) eigenvectors of (|8|) are discussed extensively in ]19|-|2l|]. It has been shown by Rashid |2l| that there is a 
one to one correspondance between each eigenvector \w) and a parent state \lm). Therefore instead of \w) it is better 
to denote a solution of (^|) as \rjlm). The relation between \lm) and \rjlm) implies a normalisation factor a; m given 
in M] and an operator such that 



7T 

\rj Im) = ai m exp(5 L z ) exp(-i- L y ) \lm) . (9) 



The parameter 8 is related to 77 by 



,x 1 ,(,)) = \/i±^, (10) 



and the eigenvalue w is expressed in terms of m by 



w = m yl V 2 ■ (11) 

Among the 21 + 1 states (||) one can identify: 

1. The states with m = which fulfill eq. (||) and (|J). These states will be taken into account and will be explicited 
thereafter. 

2. The states with m = ±1. These particular states were first introduced by Bloch and Radcliffe [Q. In a more 
convenient system of coordinates they fulfill the simpler equation with 77 = 1 [p5] , pT|| 

(L x +iL y )\w) =Q . (12) 

Therefore these states will also be considered in our paper and we will call them circular states. 

3. The states for which the parent value m is neither or ±1. These states do not coincide with the previous ones. 
However they are not orthogonal to them. Moreover they require a value of I larger or equal to 2 i.e. a tensor of 
rank at least equal to two is needed in order to generate them. We have not studied these states and it is still 
an open question to build a convenient WP by implying them. 

On the contrary the states with m = can be generated quite naturally starting from a three dimensional gaussian 
WP as shown in the Appendix and require a very simple vector operator. The states with m — have a very simple 
structure in terms of spherical harmonics that is worth to be presented shorly independently on the general solution 
found in J2~l[ . Let us denote by yL(8, 4>) these new spherical harmonics which depend on a continuous real parameter 
7? and let us expand them in terms of the usual Y m as 

y l v {6^)= cl(v)Y m (o,<P) ■ (13) 

m— — I 

The recurrence between the C l m derived from eq. (|^) implies that the sum ( p3[ ) is restricted in such a way that m and 
I have the same parity, indeed the recurrence is: 



_ ^ l{l + l)-m{m-l) 

1 - r) V 1(1 + 1) - m(m + 1) 



°m+i — °m-l 1 Z \l 77i 1 r\ Z7FZI , 1 \ ' y 1 *) 



In the following we will need the expression of y l n for I = 1 which is 
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^H (1+ ^t: a :^ d5) 



1 A3 1 



sm (/(cos (p + it] sin ( 



4 v 7T 

The combination of 9, <j> and 77 given above will be defined as 

v = sin#(cos</> + ir)sm(j)) . (16) 

One has 



W.|3$> = A ■ ( 17 ) 



Similarly, for Z = 2 the coherent states are 

3>*(M) = 
while the average of L z is 



o 9 n -l/2 

f(i + ^ 4 ) + |V 



(i+^ 2 -^(i-w + (i^m , as) 



It is interesting to point out that rj = ±1 corresponds to states with m = ±Z while the states with Z = are eigenstates 
of L x and can be the more simply written as single spherical harmonics of the angle 9' defined as 

cos 9' = sin 9 cos <j> , (20) 

and one has 

l 

y l (9 7 4>) = Yt(e', 4/) = Ct(0) Yi(9, 4>) ■ (21) 

m— — l 

The coherent spherical harmonics yt(6, 4>) have no freedom in them which allows a proper angular localization. It 
is therefore necessary to consider linear combinations. 

B. General WP. 

The most general WP solutions of (](]) which fulfill cq. (||) will be written as 

l^gcnca! = £ A< |^) . (22) 
I 

They depend on 77 which can be interpreted with the help of eq. (fl) and on weights A' which can be determined in 
order to provide a convenient angular localization. Again there will be states with 77 = ±1 that will be called circular 
and others with r\ = that will be called linear. The WP defined with other values of r\ will be called elliptic in the 



following. The justification of this name will be given in the Appendix. In section III many results will be derived for 
the particular class of WP that will be defined in the following subsection. Most of them can be seen to apply to the 
general WP ©. 

C. Exponential coherent WP. 

Among functions that are possible those which can be expanded in a power series of the variable v defined by ( |l6| ) 
are particularly interesting because they will contain all the partial waves y^. We have chosen to concentrate on the 
exponential coherent WP 
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= V 2^sinh2iV e W5infl(cOS0+i " si " ti > (23) 

which possess important properties: they fulfil eq. (g), they have direct connection to coherent states of harmonic 
oscillator (see the Appendix) and finally they have a simple geometrical interpretation. Indeed the real parameter N 
introduced there allows a proper adjustment of the angular spread. The probability density depends only on N and 
on the angle 9' defined by (|2fj|), the expression is 

If we put 77 = 1 into (18) we obtain a coherent state defined by Mostowski [^5) who wrote it as 

4>) = C~$ e N( - Sl+ia ^ a . (25) 

Here Hi and Ui are two perpendicular unit vectors (in our case U\ is along Ox and u 2 along Oy and we have eq. (|^) 
and n is a unit vector in the direction (9, cf>). 

The generalization of ( p5j ) with a parameter rj was never considered until now to our knowledge and the time 
evolution was never studied yet. Let us point out that eq. ( p3| ) can be generalized as 

%(6, <f>) = C-i e *&i+W3>a (26) 

with arbitrary but perpendicular Ui and . Calling u 3 a third unit vector perpendicular to U\ and U2 we will obtain 
WP which do not fulfill eq. (||) but rather 

ALlAL 2 2 = ±(L 3 ) 2 . (27) 

Instead of eq. (^) we would have 

(Li)=0, (i 2 )=0. (28) 

The choice of axis made in eq. ( E3|) simplifies much the interpretation and the partial wave expansion. This expansion 
will be given in subsection [II D[ 



It is not difficult to derive that our defined by (|23J) have the following average values and the following limits 
for large N: 

(%\L Z \%) = (L z ) = V [Ncoth(2N) \] ^ V [N - \] . (29) 
Therefore eq. (0) implies that 

AL y = (Ly) = i[7Vcoth(27V) - i] N — ~[N- ~] , (30) 

AL X = (Ll) = ^[Ncoth(2N) - \] "-^ £ [N - \] . (31) 

The average of L 2 Z and the total uncertainty AL 2 need independent calculations with the results: 

(L 2 z ) = {L 2 y )[l-2 V 2 }+r 1 2 N 2 , (32) 

All = <L 2 ) - (Ll) N - N I + v 2 % . (33) 
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D. Partial wave expansion of the exponential coherent state. 



We need to distinguish the case with a general value of 77 from the simpler cases with 77 = 1 or 0. (The cases with 
77 = — 1 or negative 77 are trivially deduced fom the cases with positive 77 by inversing the sense of rotation of the WP). 
For r\ = 1, i.e. circular exponential WP, or Mostowski's WP one has 



= •,(».« - ^1-^ £JSL= Y/(<>A) ■ (34) 



For 77 = 0, i.e. linear exponential WP, it is found that 



"'» ,f '"* 1 *»>"-'>-i/5SK5f W ~'' (35) 



2N 



sinh 2N 



7=0 v 



(36) 



where I I+ i(N) is a spherical Bessel function of the first kind. 

For a general value of 77, i.e. elliptic WP, the argument v can be separated into two parts 

v = sin 9 e'* + — - sin e"^ . (37) 

2 2 v ' 

Then e^" is calculated as two power series containing products of and Y l _ v . These products are then expanded in 
terms of Y^ as 

% {0, <t>) = E b lM (N, 77) l£ (0,4>), (38) 

/M 

with the weights 67 m given by 



2iV ^ (-l)'Xl + r7)]'[iV(l-77)]'' (»'00|J0)(Zn-Z'|JM) 
= V sinh(2A0 E V(20«2*')! V^TI ' (39) 



From the discussion made in subsection II A 6jm is proportional to Cjm- The known selection rules of the Clebsh- 
Gordan coefficients which appear in (39) assure that M should have the parity of /. 



III. COMPARISON WITH COHERENT STATES DEFINED IN TERMS OF BOSONS 

In this section we will compare the previous coherent states to an other set defined in terms of bosons of spin s. 
Based on Schwinger's work ]27| several angular momentum coherent states have been constructed which rely upon 
boson representation of spin s. The case with s = 1/2 was first considered by Bonifacio et al. fl33"| and studied more 
extensively by Atkins and Dobson |p6| . Mikhailov |34| has generalized this case to any spin integer or half integer. 
His work was complemented by Gulshani j35|. According to Mikhailov, a coherent state formed with 2s + 1 bosons 
depends on 2 complex numbers called a+ and a_ combined to define other complex constants otj m by 



2j 
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x jm = a^ +m aJ~ m ( , "J m ) , (40) 



for j = 0,s,2s,...,ps,... and m = -j, -j + 1,. .. ,j - 

The coherent state, called generically \as), is expressed by 



\as) = exp(-— ) Yl exp(a s/i at) |0) (41) 
= exp( — — ) E ~7=f a i m \i ms } ' ( 42 ) 

j— 0,s,...,ps,... m——j 
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|0) represents the vacuum, (/i = —s, . . . , s) is a creation operator of a boson of spin s, while the normalization 
constant depends on a + and a_ through 



n a ' = (|a+| a + |a_| a ) a '=X;W a - ( 43 ) 



The states \jms) are normalized states of angular momentum j constructed from the aj, [see Mikhailov for the full 
expressions]. The expansion ([42] ) contains integer and half integer j for s = 1/2, it reduces to integer j for s = 1, for 
s = 2 only even j occurs etc. The states (Q)-(^) are eigenstates of the annihilation operator (// = — s . . . + s) 
with eigenvalue a s ^. Mikhailov has also calculated the expectation values of various operators which are expressed 
in terms of and a M . Due to this technical simplification, the calculations of expectation values are easier than the 
work necessary to obtain formulas ( |29| ) to ( |33| ) in the case of the coherent state ^^(6, 4>). It is interesting to compare 
the properties of the state \as) to those of the general state defined by eq. ( |22] ) or to (|2^). Since a jm is generally 
non zero for all values of j and m the states \as) cannot be identified with our elliptic state for which m should 
be of the same parity as j. Clearly formulas ( |40| ) and ( |39| ) are different. At first sight one could think that (41) 
which depend on s and on two complex numbers a+ and a- describe an ensemble of WP larger than our WP (22) 
which depend only on N and rj and would therefore apply to a larger variety of physical situations. We will show 
that on the contrary the states ( pll| ) are a particular set of circular states and do not contain the freedom allowed 
with the parameter r\. Indeed if a+ and a_ are both nonzero the states (|4j) do not fulfill the conditions (|3|). In 
order to fulfill these conditions it is necessary that either a + or a_ should be zero. [In p3] it is shown indeed that 
(L x ) ~ Re(a^_ a_), (L y ) ~ Im(al £*_)]. A change of axis leads to otj m — except if m = j (if a_ = 0) or if m = —j 
(if a+ — 0). The phase of a+ or a_ can be incorporated in the phase of the states \j,m = ±j, s). Choosing for 
example a_ = and a + = k = real number the state \as) is now simply denoted as \ks) 

(k 2 ) 2s ^ 1 

\ks) = exp [— y -^— ] 22 ~n ( fc2 ) J \3 =V^m = J, s) (44) 



2 J ^ Jn\ 

j— 0,s,...,ps,... 
fc 4s 

= exp( — — ) exp[a ss at] |0) . (45) 

In the previous section, it was clear from the compact expression ( |25| ) given by Mostowski that the same physical 
state can be written by introducing three additional angles which are necessary to specify u\ and ui- The proof that 
this freedom exists also in the boson representation was given by Mikhailov [M . By a convenient choice of axis the 
state written in eq. ( |42| ) with 2s + 1 bosons and two complex parameters can be brought to the simple form ( ^ji^ ) 
with only one boson of spin s with fi = s and one parameter a ss — a 2 ^ — k 2s . Such a choice was also done by Atkins 
and Dobson p(| for the case s = 1/2. The matrix elements of the components of L given below are a particular case 
of formulas given in rcf . [p6[ . 

(ks\L z \ks) = sk 4s , (46) 
(ks\L 2 z \ks) =s 2 k 4s {k 4s + 1), (47) 

(ks\Ll\ks) = {ks\Ll\ks) = ^k is , (48) 



AL 2 X AL 2 V = -(L Z ) 2 . (49) 



These formulas and expansion (Q-[l5|) compared to formulas ( p9| ) to (|32|) show that | ks) do not coincide exactly with 
Mostowski's coherent state. For s = 1/2 Atkins and Dobson have proposed to truncate the sum over j in order to 
take into account only integer values of j. In this process the new state that we will call \k, j)j (i for integer) needs 
to be normalized properly and the formulas ( phf ) to fl48| ) cease to be valid. Eq. (ji^) still holds because the state is 
nevertheless a circular state. 

Our purpose is now to compare ^m{8, <fi) to \k, It is interesting to choose the parameters N and k in such a 
way that the limit ( p9|) for large N is valid (ry = 1). Then using ( |4^ ) one puts 



\^=N-\. (50) 
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The comparison presented in Fig. 1 shows that the probability C/ 2 to find the partial wave Y/ in and in \k, |)j 
is practically the same for N = 20. For smaller N and k (i.e. N of the order of unity) we have observed small but 
not meaningful differences. It is therefore possible, if N is large enough, to identify Mostowski's coherent state to the 
boson representation of spin 1/2 of Atkins and Dobson. 

We have shown in this section that the coherent states derived in the literature using a boson representation belongs 
to the par ticula r class of coherent circular states defined in section |n| The exponential coherent states defined in 
subsection [I C above have the advantage of a well controllable localization and also depends on an interesting new 
parameter rj. 



IV. TIME EVOLUTION OF COHERENT ANGULAR WP 

The coherent states built in section [n] are applicable to systems which have only angular coordinates on the sphere. 
This is the case of the three dimensional rotor with an axis of symmetry like an hctcronuclear molecule or some 
deformed nuclei. In the following we will assume that the eigenvalues of these systems obey rigorously the /(/ + 1) 
law and we will use the frequecy ujq written in terms of the moment of inertia J as 

= (51) 

Our WP do not allow to consider other degrees of freedom like vibrational ones or internal excitations. However 
we are confident of the interest of our work which provides a full quantum-mechanical description of the rotation of 
a pure state of a three dimensional system. 

Let us define the revival time, T rev as 

2ir 

T rev = — = (2I+l)T cl . (52) 

UJQ 



This time is twice the period of true revival of the WP. Indeed one has for the general case ( B2J ) 

{0, </>, t) general = £ h e~* ^(I+D^o tyl ^ (53) 
I 

Since /(/ + 1) is always even the period is indeed T rev /2. Nevertheless we will continue to use the same notations as 
in [JiJ. These authors have also introduced a second characteristic time T c ; called the classical time. It is defined in 
terms of the average angular momentum I defined by the average energy of the WP 

/(!+!) = £ |A/| 2 /(/ + l) (54) 
i 

For the case of an exponential WP one obtains with the help of eqs. ( f30| - |32"| ) : 

1,,, 1 " 



7(1+1) = {L?) N - {N- 1 -)+^ 



2 



N 2 - — (N — -) 
2 V 2' 



(55) 



The classical time T c \ is defined as: 



T cl = —^1 . (56) 

T c i is the period of a classical rotator having angular momentum 1 = 1. At times t = (m/n) T rev where 2m < n (m 
and n are mutually prime integers) we will use the trick developed in jlj] to write the quadratic exponential in I as 



l-l 

i 2tt I 2 JB 



« s e- l2T/ f (57) 



s=0 

It is necessary to distinguish three cases: 

a) n is odd, then I — n and all the coefficients a s are nonzero. However they have the same modulus 1/y/l. 
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b) n is even and multiple of four, then I = n/2 and the modulus of the a s have the same value as above. 

c) n is even and not multiple of four, then I — n but the a s with even s are zero, the other have their modulus 
equal to 1 / y^n/2. 

The number of values of the a s which are nonzero is called q with q = n if n is odd and q = n/2 if n is even. The phase 
of the a s can be calculated as said in Q. Using these results and inserting ( |55|) into ( |53|) one obtains eqs. (^8[),([59|) 
below: 

i-i 

9 V (6,4>, ^T rev ) general =Y A z « s e- l2 ^^ + f) y*(d,cp) (58) 

/ s=0 
2-1 

= J2^ci((>>M. (59) 

3=0 

At times £ = (m/n) T ret , any WP is a sum of g fractional WP each with a different effective time i s 

77% S 

t. = (— + j)T rBU (60) 

The fractional WP at times t s is given by 

Ki(0^,t s ) = J2 *ie- iIuot ° y^OA) ■ (61) 

There are several cases for which all the ^ s cl are clones of the initial WP defined by (^2|) for all possible valus of t s . 
There are also cases where only one of all the fractional wave is a clone for particular t s . Let us describe now these 
events keeping in mind as far as possible arbitrary A/. 



A. Cloning of circular WP 

If rj = 1 one has 

e -ii Uo t 3 yi^fi = e -i^ot s Y /ie,<f>) - Y/(9,<f>-w t.) . (62) 
Independently of the A/ the fractional waves verify the cloning property 

*2 I (0,&t.) = *(M-Wot as O). (63) 

Among all circular WP which all clone in this way, the exponential WP, which can be sharply localized in the angle 9 
by considering high enough N, clone accordingly around q directions disposed symmetrically in the Oxy plane defined 
by q values of the angle loq t s . 



B. Cloning for some particular t s 

This special situation occurs when there exist values of s such that c^o t s is a multiple of 2tt and for which a s is non 
zero. Let sq be defined by sq = n — m. One has the property 

^(0,0,i So ) = *(M,O). (64) 

This event occurs whenever n is odd or even and not multiple of four. The clone is found always identical to the 
initial WP. Obviously it is multiplied by a Sg . The existence of the clone is independent of the A/ and of rj. For the 
values of s ^ so the fractional WP are different from the initial WP. Starting from 77 = 1 their shape evolve with 77 
and we propose to call them mutants. This mutation can indeed be seen as a transition between two symmetries as 
we will show numerically in section M. 
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C. Symmetry properties of the fractional waves 



In the following discussion it will be assumed that A/ 7^ both for even and odd values of /. If some further 
symmetry is assumed (for exemple A/ = for odd I) new properties will result which will not be discussed in the 
present paper. 

The first remark is that, apart for the value sq defined above associated to a clone, the fractional waves can be 
paired for each value of m and n in the following manner: associated to s there exists an other value s' such that 

and also (noting a difference in the sign of 77 in the right side) 

^eA,tAn = Ki(0,M%. (66) 

This equation shows that the fractional waves coresponding to opposite 77, i.e. opposite (L z ), are intermixed. The 
equality ( |66"| ) is based on the equality which defines s' in terms of s 

= (mod 1) (67) 

■* rev 

from one hand and from the conjugation property 

y^ I (e,cj>) = y I _ ri *(e,<p) . (68) 

In addition to being real it is important that A/ should be an even function of 77 (like for the exponential WP defined 
by @§D). 

For 77 = the fractional waves 4^ and have the same probability density on the sphere. This leads to a 
reduction in the number of fractional waves which occur: for odd n there will be one clone plus (n — l)/2 fractional 
waves, for even n and not multiple of four there is one clone and (n — 2)/4 fractional waves, finally for n multiple of 
four there will be n/A fractional waves. 



V. NUMERICAL CALCULATIONS WITH EXPONENTIAL WP 



In this section we will describe some figures showing the time evolution of typical exponential coherent WP. The 
value of N will generally be the same and we will change the parameter 77. Most of the figures are calculated for 
N = 20. This value is typical for a rather concentrated WP. Values near unity correspond to broad WP which occupy 
the whole of the sphere and are not interes ting for our purpose. Let us remind that keeping the same N and changing 
77 produces the same probability density (24) at t = however formula (^8-39) shows that the distribution of the 
partial wave depends strongly on 77. The average I is very low for 77 = and this produce a difference in the time 
evolution which shows less structure if N — 20 and if 77 is decreased. These features can be seen if one study the 
autocorrelation function represented for three values of r) in figure |^. It is seen that this function is composed of peaks 
which have a larger width if 77 is small. The structure becomes very rich if 77 is nearer to 1. This autocorrelation 
function is very much similar to that studied in ref . p4[ . 



A. Cloning for circular WP 

The time evolution of the circular wave packet corresponding to N = 20 is shown in Fig. || and Fig. ||. In Fig. || 
a convenient set of times have been chosen to show the probability density as a function of 9 and <fi just during the 
regime of spreading of the wave packet (t < T c {) then for a few cloning times followed by the full revival for t — T rev /2. 
In Fig. H a "carpet" is shown of the section < t < T rev /2 of the probability density for 6 = n/2 which shows up to 
q = 7 clones. 

The difference between Mostowski's coherent state with that created by Atkins and Dobson \k,s = i)j is so tiny 
that it does not present any interest to be shown. Despite the analogy with the one dimensional results of the infinite 
square well [Q and the two dimensional rotor jll| some new interesting aspects of our results need to be stressed. 
Indeed the circular WP spread in a way in the <f> direction and clone around the Oxy plane, which is natural since 
there is a linear momentum along Oy initially. However there is no change in time in the azimuthal spread. 

The cloning mechanism found in quantum mechanics is not possible for a single classical particle however it will 
appear if one uses the ensemble interpretation of quantum mechanics as underlined by authors of ref. B . 
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B. Linear WP 



For 77 = the fractional WP has cylindrical symmetry around Ox at all times since it is written as 

¥ d = E 6, e - 2 " /t/T "" Yj{e',4>'). (69) 

The quantity 27rsin6*'|^' c /| 2 is represented in Fig. |[ This shows that the revival wave packets are rings on the 
sphere and, due to this special topology do not clone the initial wave function. 

The time evolution of SS>o(0' ,4>' ,t) is represented in Fig. |[ Since the wave packet is in fact one dimensional, the 
"carpet" representation provides the essential features of the time evolution. However in order to produce a similar 
richness than for 77 = 1 (and a similar expectation value of L 2 ) we have increased the value of N for 77 = to N — 50. 
On a sphere areas of constant probability density could be represented by parallel circles centered on Ox axis. The 
pattern of the carpet shown in Fig. || is very much resemblant to that discussed recently for a special wave packet in 
one dimensional box in ref . [B8fl . Indeed there is a superposition of ridges and valleys with simple slopes as a function 



of t. The interpretation of this effect can be given in similar terms as in 38 1. Note that the quantity plotted in Figs.|6j 
and H has the same boundary conditions as the wave packet on the edge of the box. This produces a reflection effect 
in Fig. H totally absent in Fig. [| since the boundary conditions are different on the circle. Thus for 77 = the WP 
spreads uniformly in all directions defined by the angle <j)' . 

Such a WP is full nonsense for a single classical particle and takes sense only with the ensemble interpretation. 

Another interesting linear WP corresponding to 77 — * 00 but keeping 77 TV finite is defined as follows: 

V47T 



It is derived from the harmonic oscillator WP of the appendix by keeping only the term in p y . This WP has its 
probability density uniformly distributed over the sphere and obeys the equation 

L y y , qN = . (71) 
It has therefore cylindrical symmetry around Oy and depends on the angle 6" defined by 

cos6»" = sin6» sin./). (72) 
Its expansion in spherical harmonics contains now spherical Bessel functions 

9 vN = V^T+Tjj^N) jtf (0" , 0" ) . (73) 

Obviously the fractional waves have also cylindrical symmetry around Oy but there is in addition, for t = t So , a 
uniform clone which interferes with all the other fractional waves. The time evolution is shown in Fig. ^ for rjN = 20. 
For the small values of t, like t — 1/100 T rev , the WP is concentrated almost totally on the hemisphere with y > 
with a spike along Oy surrounded by concentric rings. For t = l/50T rev the same behaviour occurs but this time on 
the hemisphere with y < 0. For other times both parts of the sphere are covered with rings and the spike also occurs 
on both sides of the Oy axis. For t = T rev m/n with small m/n a symmetry between the two hemisphere takes place. 
The existence of the clone can be seen clearly as a small uniform background at times t = l/25T rev , t = l/10T rev . 
There is always a strong interference between the fractional waves which does not allow to make a clear counting even 
for small values of m/n. 

C. Elliptic WP 

For the general elliptic wave packet, as deduced from the previous discussion there are no clones, but partial revivals 
with different topology. Due to the change of symmetry, it is indeed necessary to make a smooth transition between a 
system of clones located for rj = 1 in the Oxy plane and a system of rings discussed in the previous subsection. In the 



system of coordinates adopted in subsection II A and corresponding to 77 = these rings have Ox as symmetry axis. 
The transition from the clones for 77 = 1 to the rings for 77 = is made by developing for 77 smaller than 1 a system 
of pairs of crescents perpendicular to the Oxy plane. This transition is clearly visible for particular fractional revival 
times in Fig. H. For a very small value of 77 the upper part and the lower part of one crescent meet the corresponding 
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parts of a symmetric crescent in order to build up such a ring symmetric around Ox. Again this change of topology 
of this construction forbids the cloning of all fractional waves. Most of them revive in shapes different from that of 
the initial wave packet. We propose to call these fractional wave packets mutants. Example of time evolution of wave 
packet with 77 = 0.5 is given in Fig. |. It is well seen both in Fig. [| and in Fig. |^ that clones and mutants can occur 
in the same time. For example for f] = j an d t — | T rev there is a fractional revival at 9 = ^ and (j) — with the 
same shape as the initial wave packet and two crescents which almost close. This situation is similar to that which 
exists for 77 = also for t — ^ T rev . For T] — \ and t — ^ T rev the two crescents do not form a ring. Among the 
three fractional waves which are present there are two which are identical to each other but with larger spread in 
9. For this value rj = i, there are three different topologies for the five fractional waves. In other publications on a 
different system pOfl , we had already found the transition from a gaussian three dimensional wave packet to a vortex 
ring. Refs. |50j were devoted, as well as our previous works, for the time evolution of such coherent waves in the case 
where the hamiltonian contains a spin orbit interaction in addition to the harmonic oscillator potential. If the spin 
is oriented along the initial wave packet displacement (Ox axis), the cylindrical symmetry is imposed on the system 
and preserved during evolution, a vortex rings appear. 

The transition between the two shapes of fractional wave packets can be explained by comparing, as a function of 
77, the spread of the angle 9 to the spread of the angle <p. If 77 = the wave packet has no privileged direction on the 
sphere. It spreads therefore equally. In the case when r/ = 1 the wave packet is peaked near the plane with 9 = j 
and the spread in 9 is strongly reduced. The spread in </> is seen in the scenario of cloning in the plane xOy. For 
intermediate values of 77, there is a competition between the two effects which manifests itself most strongly in the 
shape of the fractional waves. 

For 77)1 (results not presented) we have observed a strong reduction of the spread in 9. For particular values of 
(m/n), some mutants can be peaked with higher amplitude than their neighbours. This is in direct connection to the 
increase in po which was necessary in the initial wave packet. 



D. Final remarks 



We can measure the aperture of the probability density of the WP by the solid angle SI = 47r/(47V+ 1) corresponding 
to the cone defined by tan(#'/2) = l/2\/N. The time of spreading t v of this wave packet of width AL V ( |33| ) is of the 
order of 

t. (74) 

The difference in spreading with 77 is clearly seen in the autocorrelation functions represented in Fig. |[ We can also 
define maximum number of fractional wave packets that can be observed knowing that tan(6*'/2) = 1/2^/N as 

(75) 



arctan(l/2vA0 

This number is confirmed by the observation of the 'carpet' of Fig. |J for N — 20 for which we have up to 7 clones. 
The lifetime r' \ of a system of q fractional wave packets can be estimated as 

coo t'„ AL v = — . (76) 

q 

Also in Fig. |i| one sees clearly for q = 2, 3 the large intervals of time during which clones are observed. 



VI. COHERENT STATES FOR A SYMMETRIC TOP 



It is natural to enlarge our previous study to rotational coherent states of symmetric tops which contain an additional 
degree of freedom. There are several possibilities in the literature for constructing such a state. Our aim is again 
to choose the axis of coordinates in order to eliminate irrelevant parameters and to make evolution understandable. 
This is possible only in the case of symmetric top. We will assume that the moments of inertia in the intrinsic frame 
Jx = Jy and introduce the parameter 5 such that: 
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The energy spectrum is then: 



E(I,K) = £-[I(I + 1) + SK 2 }, (78) 

A convenient rotational wave packet was defined and studied by Janssen [f28| on the basis of the work of Perelomov 
36| and Schwinger p7[ . This wave packet denoted by three complex numbers x, y and z is a mixture of D I MK functions 

\xyz) = exp(--yy*(l + xx*)(l + zz*)) (79) 




(2/)! x I+M y 2I z I+K \IMK) 



{I + M)\{I - M)\{I + K)\{I - K)\ 



The sum over I contains integer as well as half integer values of I i.e. this state extends that presented by Atkins and 
Dobson. On a similar line as in section III we call \xyz)i its normalized projection onto the space with integer /. A 
particularly convenient choice of x, y and z as well as of the system of coordinates decreases the number of parameters 
to two and reduces the summation to I and K only. The simpler wave packet is then: 

fain hV+ K (cn<i ^V~ K 

\r, A) - exp(-r) Y (-1) I+K (2rf [ 2> |z _ m) (80) 

If Lx, Ly , Lz are the components of L in the intrinsic axes, L Xl L y ,L z those related to the laboratory axes, one 
has according to Janssen: 

(r\\L x \r\) = (r\\L y \r\) = 0, (rX\L z \rX) = -r , (81) 



(rX\L z \rX) = -rcosA, (rX\L x \rX) = -rsinA, (rX\L Y \rX) =0, (82) 

(rX\L 2 \rX)=r(r + ^). (83) 

Since the WP given by ( |80| ) contains only components with M = —I it is a circular WP with r) = — 1 which fulfills 
cq. (||) with the components of L taken in the laboratory system. From Janssen's work the components in the rigid 
body system verify the equation 



while their product takes the value 



AL 2 X = ALy , (84) 



AL 2 X AL Y = . (85) 

4 cos^ A 



Therefore eq. (0) is also verified for the components in the rigid body system if OZ is directed along L. On the 



other hand the expansion (80) of the rotational coherent states in terms of \I — IK) has the same coefficients as the 
expansion of Atkins-Dobson in terms of angular momentum eigenstates \IK). The two expansions are connected by 
defining a+ and a_ as: 

a + = y/2r sin ^ , a_ = y/2r cos ^ . (86) 



The projection |r, X)i deduced by restricting (bfj) to integer values of I and K verify eqs. (|8lj) to (|S3j) only 
approximately. However if r is large enough, these equations are obtained in a very good approximation. The time 
evolution of the wave packet |r, X)i will be studied in the next section for a rigid body symmetric rotor. 
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VII. TIME EVOLUTION OF JANSSEN'S COHERENT STATE 



The energy spectrum of the axially symmetric rigid rotor is written as: 

E IK = huj [1(1 + 1) + SK 2 } . (87) 

We will study the time evolution of a wave packet deduced from ( pp[ ) with uj = h/(2J x ) and with average values 

(rX\L z \rX) = K ~ -rcosA (88) 
(r\\L z \r\) = I ~ -r . (89) 

The energy Eik is written by taking K and I as reference, and defining k\ and k 2 as: 

kx=I-I k 2 =K -K . (90) 

E klk . 2 = 7kj [I (J + 1) + 5K 2 } + hu (21 + l)ki + hu 5(2K)k 2 + huj Q k\ + Huj q S k\ . (91) 

We now follow the lines drawn by Bluhm, Kostelecky and Tudose |[o| who have considered the time evolution of a 
system which depends quadratically on two quantum numbers, in our case I and K . There are 4 time constants; the 
first pair defined as 

T '' = Ithvi T - = - = (2/ + 1) 1^ (92) 

is related to the motion around Oz axis (laboratory axis), i.e. it is connected to the Euler angle a, the second pair 
plays a parallel role, it concerns the motion around the symmetric OZ axis and the Euler angle 7 

K 2tt 1 2 1 + 1 j K 2tt 1 r 

cl uj S2K 6 2K cl rev Suj S rev ' { ' 



The system of coordinates and the parametrization in (80) enables to profit by the separation of variables in the 
state \I — IK) since: 

(a(3j\I - IK) = e laI dL IK (P) e^ K = Dl IK (a, 0, 7) . (94) 



The wave packet QSOj) at time t with conditions (|8q ) and (g9J) and integer / and K will be denoted as \I K) and 
written as: 

(aHlK)t - ]T C IK (r,X)dL IK (0)e i ^ I - 2 ^ I + 1 ^/ T r^ e -KiK+2.K- t/T«j (Q5) 

IK 

= e-^^ £ C klk2 (r, A) dU W (96) 

x ^[a^^ikt/T^+kf/T^t] e -i[-rk 3 +2n(k 3 /T%+k%/T? ev )t] _ 

The summation on / and K has been changed to a sum over k\ and k 2 and the coefficient Cjk &-ik nas been given 
the new indexes. The discussion of the time evolution of ( p5[j9^ ) follows in a straightforward manner that given by 
Bluhm, Kostelecky and Tudose frL0[. The crucial parameter is 5. If T^ ev and T^ ev are not commensurate there is no 
cloning, however for t = (m/n)Tj^ there are partial revivals in the variable a: i.e. the wave packet is a superposition 
of q fractional wave packets peaked regularly along the Oz axis, (q = (n/2) if n is even, n in other cases). For 
t = (m'/n')T£ v the same scenario of partial revival occurs but this time there are q' fractional wave packets peaked 
around the OZ axis, [q' = (n'/2) if n' is even, n' otherwise]. 

The interesting situation of commensurability of T^ ev and T^ ev allows, on the contrary, the construction of q 2 clones 
for all the time such that 

m 1 K 

t T rev (97) 
n 

The revival time is the least common multiple of T^ ev and T~^ v 

vT 1 =rT K — T I,K (98) 

F^rev rev rev ' x^^J 
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The time evolution of the rotational wave packets is presented in Fig. |l^ and Fig. 111]. 

The particular choice K — leads to Tj" = oo but T^ v is finite as well as and Tl ev . For the smaller value of 
i, the behaviour of the wave packet around Oz and OZ is different as seen from Fig. O. There is indeed a classical 
rotation and spreading around Oz while no rotation occurs around OZ, only the spreading is observed around this 
axis. 

Fig. [l^ illustrates the case of an irrational value of 5 where there are partial revivals for times t = (m/n)T^ ev 
around Oz and for times t — (m! /n') T^ ev where the revivals are around OZ . Note that the concentration of the wave 
packet for some definite values of a (or 7) has no influence on the other variable 7 (or a) respectively. 

Fig. [ll] illustrates the case of a rational value of <5, for which there are q 2 clones. The proof that there are q 2 clones 
if condition (^8|) holds is rather simple, one applies twice the method of Averbukh and Perelman to linearize the 
exponential containing k\ and k 2 in eq. (p5[-p6[). 

Finally we want to stress that the variable (3 does not play a role in the time evolution. Obviously this is due to 
our particular choice of axis in eq. (80). The role played by (3 for the axial rotor is similar to the role played by 6 in 
the diatomic molecule. 



VIII. CONCLUSIONS 



We have discussed the time evolution of angular momentum coherent states that are built from intelligent spin 
states, but as pointed in the Appendix, can be constructed by a simple, purely geometric, generating procedure. The 
basic ingredients arc the three dimensional gaussian wave packets of the harmonic oscillator. We have used two of the 
parameters of the gaussian to derive an ensemble of coherents states which depend finally on two parameters N and 
rj. The angular distribution of the probability density depends only on A, while the momentum distribution depends 
both on A and 77. The uncertainty relation (j|) is valid for all A and 77. Such a variety of wave packets do not arise 
from the previous works on angular momentum coherent states with bosons. It is for example simple to show |32j ] 
that there is no place in the Atkins and Dobson or Mihailov's work for the wave packets (^0-41) which correspond to 
77 = 0. Also the coefficient a.j m defined by (^) in terms of a + and a_ cannot be put generally in the form such that 
( [l4|) which ensures that biM is zero if / and M have an opposite parity. 

In the framework of the scenario of ref. [Q we have shown that for 77 = 1 the wave packets spread at t = 21 T rev into 
clones. However, the difference between the symmetries for 7/ = and r\ = 1 necessitates a change in the topology of 
the fractional wave packets. This change is such that these wave packets can be named mutants after a certain amount 
of deformation of their shape. Aronstein and Stroud Q have shown "that, the infinite square well is an ideal system 
for fractional revivals since all wave functions exhibit fractional revivals of all orders" . The same statement applies 
to the case of the diatomic molecule but is limited however to circular wave packets, i.e. to those which contain only 
M = I partial waves in a convenient system of axis. The generation mechanism used in this paper can be extended 
to shapes different from gaussian distribution. Clones will always result in this particular condition. For the wave 
packets which contain populations of sublevels M ^ I, mutants will appear. All these results hold independently on 
the realization of the uncertainty condition (^). For systems where energy is quadratic in one quantum number we 
have shown that the time evolution can exhibit a rich structure of cloning and mutation of fractional revivals. It is 
a challenging question to find a proper way to excite such wave packets in a rotational band of a molecule or of a 
deformed nucleus. Let us briefly describe the efforts made under this perspective. 

The wave packets with 77 = or wave packets with M — have been already considered several times. In a reference 
work on Coulomb excitation of nuclei with heavy ions R.A. Broglia and A. Winther |3l[] have shown that such wave 
packets are generated in backward scattering. Their population of the excited states of a rotational band with K = 
do not coincide with the very simplified expression ( 40|^4l| ). Indeed interference effects during the excitation process 
modulates this population to a large extent (Fig. 8, p. 82 of ref. plfl). However, we have studied the time evolution of 
such a wave packet and this will be the subject of a future publication 40 1. 

The authors of ref. p2| have written an important analysis on coherent rotational states, completing the work 
of [Q. They have already used the wave \k, |)j of section III derived from [ p6| |, and they have also considered an 
extension toward other symmetries like Sp(2,R), Sp(4,R). Time evolution was studied but only for the average values 
of operators like quadrupole moments. The physical systems considered were the nucleus 238 U and the molecule CS2. 
The time evolution of expectation values of more general operators for wave packets due to quadratic dependence 
was studied in p7| ] . The operators studied in ref. p2| ] have very strict selection rules and the time evolution reflected 
there is much poorer than what is exhibited by the full wave packet and by the average values of the more complex 
observables |57| related to the position. It is, however, clear that all the wave packets considered in ref. |32| must 
evolve according to the scenario of ref. jjj with a very rich sequence of changes of shapes. 

Our work comes very near from the spirit of two recent publications devoted to molecular physics. In ref. |39t] the 
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population of a set of rotational states of a molecule by an intense laser has been calculated. However, the number 
of states excited in this way is rather small, I < 4, the symmetry implies only states with M = and therefore the 
WP produced is not concentrated on the sphere as sharply as ours with TV = 20 or 50. In ref. O] Ortigoso has 
studied how to tailor microwaves pulses in order to create rotational coherent states for an asymmetric-top molecule, 
which are built from Radcliffe's intelligent spin states discussed in section It is gratifying that such a mechanism 
which involves optimal control theory is possible. However our demand is more ambitious since one needs to combine 
intelligent spin states of different angular momentum in order to achieve a proper angular concentration. In addition 
one has to find a mechanism that allows to change continuously our variable rj, i.e. the relative uncertaities in L x and 
L y , in order to explore the whole variety of our states as well as of the intelligent spin states. It is indeed interesting 
that Radcliffe's states have been used quite thorougly for example in re f. I|42| ,[43| but our paper points toward the fact 
that a richer structure exists nearby in accordance to older works p9|-pl|. As said in the text we have not used the 
larger class of intelligent spin states in the case of the diatomic molecule and we are also conscious that our work on 
the symmetric top leaves open possibilities of coherent states that have not yet been explored. 

It is worth to mention very recent paper by Chen and Yeazell 0] on an analytical wave packet design scheme that 
is able to create desired Rydberg wave packets and control their dynamics. 

We believe that we have enriched anyway the exemples given in the reference ]To| by pointing out how different the 
fractional revivals may be. 



APPENDIX A: CONNECTION BETWEEN THE EXPONENTIAL COHERENT WP AND COHERENT 

STATES OF THE HARMONIC OSCILLATOR 



The exponential WP defined by eq. ( |23| ) can be manufactured from gaussian wave packets in three dimensions 
which are coherent states of the harmonic oscillator. Such a general WP with an average position ro and momentum 
po are written as 



1 



l0 ,i i exp 



(f-ro) 2 .pb-r 
+ I 



2er 2 ' " h 

One gets rid of three unecessary parameters if one chooses the axis in such a way that 

p 0z =0, r = xr . 

With such choices the probability density is 



1 



■ exp 



2r r 

exp | — — sin 6/ cosi 



(Al) 



(A2) 



(A3) 



(2tt) 3 / 2 ct 3 

Apart from a normalization factor it coincides with eq. ( p4| ) if r and ro are both chosen as 

r = r =(J VN . (A4) 

The solid angle f2 = Air /(AN + 1) is thus the angle with whi ch the width a of the density of the gaussian WP is 
observed from the center of the sphere of radius given by ( A4). This choice still leaves two parameters po x and po y 
free. The choice which leads to (E3ft is 



Po,, 



0. 



roPo, 



ri N h 



If one takes the width of the harmonic oscillator a = \ff\Jmu) the value of po y given in (A5) leads for rj = 1 to 



p 0y — m tq uj 



(A5) 



(A6) 



The coherent state (Al) associated to this value evolves around a circular trajectory in the field of the harmonic 
oscillator. In the same manner the value rj = is associated to a linear trajectory and the other values of ij correspond 
to elliptic trajectories. In the second case the initial point is either the apogee or the perigee of the ellipsis. In this 
manner more general WP can also be constructed which start from arbitrary points according to a nonzero value 
given to p x0 . Let us define a parameter e by 



r Po x = eNH . 



(A7) 
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The exponential WP fl23Q becomes in these new conditions 



VT/ (/-;,,) Afsinflf(l+»e)cos0+^sin0| \ s 

w ^^'^-\/ 27rs inh2iV 1 ' 



and verify the condition 



[(1 + ie)L x + irjLy] = . (A9) 



The probability density is again given by ( p4| ) . However these more general states do not fulfill condition (Q) . Indeed 
for those states which obey this equation there exists a value of r\ given by (0) such that eq. (|^) is obeyed and eq. 



(AE) is not. It was then consistant to consider only those WP with e = 0. This result is in accordance with an older 
result by Rashid plj and mark a difference between so called quasi-intelligent spin states which solve (A8) but not 
satisfy (H) and the intelligent spin states for which e = 0. 
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FIG. 1. Probabilities to find the partial waves Yf in the coherent states of Mostowski (solid line) and Atkins-Dobson (dashed 
impulses) for parameters N = 20 = l/2(k 2 + I) ensuring the same angular velocity for both wave packets. 

FIG. 2. The autocorrelation function for N = 20 and different values of parameter rj corresponding to a smooth transition 
between two different symmetries. 

FIG. 3. Time evolution of Mostowski's wave packet with N = 20. The left column represents changes of the probability 
density during short term evolution, the right one at fractional revival times. The probability density at times equal 1/5, 1/3 
and l/2*T rev are identical with those presented at 1/10, 1/6 and n*T rev , respectively. (The vertical scale is not the same in 
all figures.) 

FIG. 4. Time evolution of Mostowski's wave packet with N = 20. The \^/ m{9, 4>,t)\ 2 f° r fixed 6 = ty/2 is presented in the 
contour plot. The larger values of \^>m(6, cf>, t)\ 2 result in heavy lines due to almost overlaping cuts for fractional revivals. One 
can clearly see the fractional revivals of orders 1/7, 1/12, 1/5, 1/8, 1/6, 1/4 and so on (corresponding times are 1/7, 1/6, 1/5, 
1/4, 1/3 and 1/2 of Trev)- 

FIG. 5. Time evolution of the ty c i ( |69| ) for N = 50. The probability density 2tt sin | \& c ; | 2 is presented as a function of 6' 
and t within one revival period. 

FIG. 6. Time evolution of the wave packet (|35]-p6[) for N = 50. The probability density 27rsin#'|*I'| 2 is presented in the 
contour plot. The larger values of 2-k sin 0'|\E'| 2 result in heavy lines due to almost overlaping cuts for fractional revivals. 

FIG. 7. Time evolution of the WP @ with infinite r) but t)N = 20. At t = the probability density is 1/4tt. See 
explanations in the text. 

FIG. 8. Transition of fractional wave packets from exact clones (r/ = 1) through developing crescents (rj — 1/2, rj — 1/4) to 
ring topology (j) = 0) is demonstrated for two fractional revival times t = 1/3 * T rev (left) and t = 1/4 * T rev (right). The 
fractional waves called mutants are clearly seen in the lower rows of the figure. 

FIG. 9. Time evolution of the coherent state deduced from elliptic motion (rj = 0.5). The time sequence presented here is 
the same as in Fig. []. Comparing to this figure the crescents are clearly visible, they are the most pronounced for t = j T rev . 

FIG. 10. Time evolution of Janssen's coherent state for axially symmetric top ( p5|j9^ ) for irrational 6 = l/\/3 implying 
T/ e „ = 1/V3T^, V ~ 0.577T/L- Shown is the probability density \{af3j\IK)\ 2 for (3 = tt/2, 7 = 4 and K = 0. Notice that the 
vertical scale is not the same in all figures. 

FIG. 11. The same as in Fig. [l^ but for rational 5 = 1/2 implying = T, K ev = 2T/ e „. Clones for times 1/6, 1/3 and 

1/2 Tfew are clearly visible as well as full revivals. Notice that the vertical scale is not the same in all figures. 
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FIG. 1. Probabilities to find the partial waves Y/ in the coherent states of Mostowski (solid line) and Atkins-Dobson (dashed 
impulses) for parameters N = 20 = l/2(k 2 + 1) ensuring the same angular velocity for both wave packets. 
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FIG. 2. The autocorrelation function for N = 20 and different values of parameter r\ corresponding to a smooth transition 
between two different symmetries. 
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FIG. 5. Time evolution of the 9 c i (69) for N = 50. The probability density 27rsin#'|^' c ;| 2 is presented as a function of 8' 
and t within one revival period. 
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